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Abstract 

We look back at early efforts to approximate the large Nc Feynman diagrams of QCD as 
very large fishnet diagrams. We consider more carefully the uniqueness of rules for discretizing 
P"*" and ix'^ which fix the fishnet model in the strong 't Hooft coupling limit, and we offer 
some refinements that allow more of the crucial QCD interactions to be retained in the fishnet 
approximation. This new discretization has a better chance to lead to a physically sensible 
"bare QCD string" model. Not surprisingly the resulting fishnet diagrams are both richer in 
structure and harder to evaluate than those considered in older work. As warm-ups we analyze 
arbitrarily large fishnets of a paradigm scalar cubic theory and very small fishnets of QCD. 



Copyright 2000 by The American Physical Society 



^E-mail address; beringSphys .uf 1 . edu 
E-mail address; rozowsky@phys.ufl.edu 
^E-mail address; thorn@phys.ufl.edu 



1 Introduction 



With all the recent effort devoted to the search for a solution of large Nc QCD Q as a classical 
string theory Q, it is appropriate to reassess earlier efforts to accomplish this goal. In this article 
we wish to refine and extend the formulation and calculational methods developed in the effort of 
the late seventies to systematize a fishnet Q approximation to large Nc QCD @, ^. 

The larger goal here is to set up a discrete model of infinite Nc QCD which, when analyzed in a 
weak coupling expansion [NcQ^ <^ 1), reproduces perturbative QCD and, when analyzed in a strong 
coupling limit {NcQ^ — > oo), describes what we choose to call a "bare QCD string". Since QCD is 
supposed to confine at all values of the 't Hooft coupling, the infinite Nc glueball should actually 
be noninteracting over the whole range of couplings. However, its composite internal structure is 
generally expected to be quite complicated, and it is only in the strong 't Hooft coupling limit that 
the internal structure of an infinite Nc glueball can be as simple as that of the "fundamental string" 
of string theory. We regard it as an open question whether the bare QCD string can be identified 
with one of the known fundamental strings or is an entirely novel object. We hope that our efforts 
will eventually settle this issue. We tentatively identify the bare QCD string with the object whose 
propagation is described by the so-called fishnet diagrams. 

As shown in Q the fishnet diagrams by no means exhaust the planar diagrams of 't Hooft's 
Nc ^ oo limit. Fishnets are certainly planar, but they are also very large in both directions: 
there are many lines and many interaction vertices. It is natural to try to associate such diagrams 
with strong 't Hooft coupling Ncg'^ oo, but as with all strong coupling expansions one must 
first define a cutoff theory which controls the size of the kinetic energy of the system. In Q the 
choice was to evaluate all graphs on a light front and simultaneously discretize = (P'^ + P^)/\/2 
and x'^ = {x^ + x^)/^/2. In this model the large 't Hooft coupling limit singled out large fishnet 
diagrams whose continuum limit was a seamless world sheet. As usual with strong coupling limits 
this conclusion is highly sensitive to the cutoff model. 

What one hopes when resorting to strong coupling methods is that although the limit strongly 
distorts the quantitative details of the continuum theory, the qualitative physics is shared between 
the continuum and lattice models for all couplings. In standard lattice gauge theories this hope is 
usually expressed as requiring that the lattice model exhibit no phase transition as the coupling 
constant is varied from strong to weak coupling. Probably the most familiar case in which there is 
a phase transition is "compact QED" whose strong coupling limit shows confinement, but whose 
continuum limit is a theory of free photons. 

Although the existence of a phase transition at finite coupling is usually extremely difficult 
to detect, it is the case that in some situations our lattice fishnets can be seen to be completely 
irrelevant to the physics of large but finite coupling. In |Q] this possibility was noted in the context 
of scalar Xcj)'^ theory. The qualitative physics of the seamless fishnet diagrams is that the quanta of 
the field theory are bound into a linear polymeric chain. However, one can examine at next order 
in the strong coupling expansion the nature of the interaction that should be responsible for this 
binding. For A > this interaction is repulsive, and the seamless world sheet given by the strong 
coupling limit is a purely formal artifact. In contrast, for A < the interaction is attractive, and it 
is qualitatively correct to imagine that the nearest neighbor quanta form very tight bonds in the 
strong coupling limit. 

A serious shortcoming of the QCD fishnet model attempted in ||5| is that the basic gluon-gluon 
quartic interaction retained in the strong coupling limit favored the alignment of the gluon spins. 
This defect was not apparent, however, because the leading fishnet structure was explicitly an even 
function of this interaction and in fact described an anti-ferromagnetic spin arrangement. The 
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inherent instability of the system would only be seen at non-leading order. We aim to improve this 
situation in the present work by proposing a discretized model whose fishnet approximation retains 
both the "contact" interactions, with their ferromagnetic tendency, and the one gluon exchange 
interactions with their anti-ferromagnetic tendency. 

The rest of the article is organized as follows. In Section^ we give a self contained review of the 
Feynman rules in light-cone gauge as well as the discretization rules largely following ^ . However, 
we treat the longitudinal modes differently. We represent the "induced quartic interaction" of light- 
cone gauge by the exchange of a short-lived fictitious spin quantum: its propagation is limited to 
a small number of discrete time steps. This idea motivated another departure from Namely, we 
also choose to represent the basic quartic gluon interaction by the exchange of another short-lived 
fictitious spin quantum. In this way all vertices of the discretized Feynman rules are cubic, and 
are accordingly all treated on the same footing in the strong coupling limit. In Section |3| we show 
our discretization in action by computing the gluon self energy at one loop order. We see how 
the ambiguities inherent in spreading out the quartic vertices begin to be resolved by requiring 
Lorentz invariance. The propagators of the fictitious scalars are multiplied by fk, h^, where k is 
the number of time steps, and = /i^ = 1. Lorentz invariance of the self-energy constrains 

the moments X^^/fc/^ and X^fc/ifc/^- In Section ^ we describe the fishnet approximation. As a 
warm-up, we give a complete analysis of the leading fishnet diagrams of a paradigm matrix scalar 
field theory gTrc/)'^. Then we describe the more complicated situation of QCD. We do not attempt to 
analyze the arbitrary QCD fishnets here. Instead in Section^ we study the sum of planar diagrams 
for small values of M, the number of units of carried by the evolved system. Finally in Section 
^ we collect some concluding remarks and sketch future directions for this program of research. 



2 Light-cone Feynman Rules and Discretization 
2.1 Propagators 

The gluon propagator in light-cone gauge A- = is given in momentum space by 
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Evaluating the p integral leads to the following expressions for the individual components of D^'^: 
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where the arrows indicate the imaginary time versions (r = ix^). In this paper latin indices will 
always refer to the transverse components. We shall also find it convenient to use a complex basis 
for transverse indices, defining = {V^ + iV^)/V2 and = {V^ — iV'^)/\/2. In this basis the 
metric has values r/vv = i]aa = and ryvA = ^av = 1- 

Here we are assuming that A~ = —A^ has not been eliminated from the formalism. Since Gauss' 
law relates to the transverse components through a constraint not involving time derivatives, it 
is possible to explicitly integrate A^ out (see for example [P) leaving the transverse components as 
the only independent variables. In that case the Feynman rules would only employ the transverse 
propagator D^'K Graphically, one achieves the same result by showing that the contributions of 
and D lead to modified cubic vertices and a new induced quartic vertex which arises from 
the 5{x^) term in D . 

2.2 Vertices 
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Figure 1: Cubic and quartic gauge vertices for imaginary time. 



We shall present the primitive cubic and quartic vertices as 't Hooft did in his presentation of 
the 1/Nc expansion ||^. Since the vertex assignments lack permutation symmetry, it is understood 
that all permutations of them must be included. The double line notation makes clear what powers 
of Nc must be included with each topology. However, since we shall be dealing exclusively with the 
planar diagrams of the Nc ^ oo limit, we shall dispense with this refinement in order to reduce 
clutter. To correctly use these rules at finite Nc, the double line notation should be restored. With 
this understanding, the primitive cubic vertices are given by 



r^^^ = r^^^ = 

= -i9{Q2-Qir- 

= -i9{Q2-Qir- 

= r^''^ = +^9{Q2 



and the primitive quartic vertices by 



tAAVV 
4 

tAVAV 
4 



-9{Q2-Qif 
-9{Q2-QiY 

Qi)^ ^+9{Q2-Qiy, (4) 



= -2ig'^-2g\ (5) 



3 



where the arrows indicate the appropriate vertices to use with imaginary time. In hght-cone gauge, 
only transverse gluons participate in the quartic vertices. Our convention will be that all momenta 
flow into the vertex. Also, the index V (A) will be represented graphically by attaching an outgoing 
(incoming) arrow to a line (see Fig. |^). The ordering of indices will be counterclockwise around the 
vertex. As usual, each vertex is associated with an integration over x"*" and conserves the transverse 
and + components of momentum. Then each unconstrained momentum is integrated with measure 
(fpdp+/(2T:f. 



2.3 Discretization 

To give a nonperturbative model for the summation over planar diagrams, it was proposed in Ref. 
1^] to simultaneously discretize p"*" = Im and imaginary time r = ix'^ = ka, with k, I running over 
all positive integers. The use of imaginary time converts all oscillating exponentials to damped 
ones, and removes all i's from the Feynman rules. Thus the i occurring in each vertex is combined 
with the dx~^ to form dr. (Because of time translational invariance, the time integral for one vertex 
in each connected diagram should be omitted, leaving one factor of i unabsorbed. Conventionally, 
we shall omit this last i in the quantities we calculate.) 

It can be easily seen that the lattice constants a, m only enter the sum of graphs in the ratio 
Tq = m/a. First notice that only this ratio appears in the exponents. Since each propagator is 
nominally integrated over its p"^, J dp'^ provides a factor of m to cancel one from a l/p"*" 

prefactor in each propagator. Further J — > a ^ at each vertex which also has a nominal p"*" 
conserving delta function that supplies a 1/m, so each vertex supplies a factor I/Tq. Finally, every 
— index of a propagator will be matched with a + index of a vertex, which will involve a factor 
of p+ and hence supply a factor of m to cancel the extra factor of 1/m in D^^ and to convert the 
extra factor of 1/ma in D to Tq. 

The discretization of D involves some ambiguity in the interpretation of the term involving 
6{x^). With x~^ continuous, this term collapses the two cubic vertices it connects into an instan- 
taneous quartic interaction local in x but p^ dependent and hence nonlocal in x~ . Indeed this is 
precisely the well-known quartic vertex induced by elimination of Aj^ in the Hamiltonian formu- 
lation of light-cone gauge. In this approach the remaining part of D combines nicely with the 
contributions of to yield a modified cubic vertex for transverse gluons only: 

f = -25 f {QtQ^2 - QtQXl (6) 



Qi Q2 

which are the vertices appropriate to imaginary time. With longitudinal gauge fields completely 
eliminated in this way, discretization could then proceed as usual by discretizing the x^ and p~^ 
parameters of the transverse gluon propagator. In addition, one has to exclude, in some ad hoc 
manner, the p^ = exchange part of the induced quartic vertex which is infinite as it stands. 



Moreover, the set of "tadpole diagrams" necessarily excluded in our discretization (see Sec. 2.4 
below) is enlarged by the induced quartic interaction and, since the new quartic interactions depend 
non-trivially on p~^ , the p"*" dependence of the necessary counter-terms will be more complicated. 
Nevertheless, an attractive feature of such a treatment is that the modified cubic vertex is manifestly 
invariant under the light-cone Galilei group: Q — > Q -|- Q^V. 

We shall follow a different path, more in the spirit of the sum over histories. The idea is to 
exploit our discretization of x"*" to give a more fiexible interpretation of 6{x^), which retains a 
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Gaussian damping factor and maintains Galilei invariance throughout. This can be done by the 
replacement: 



A:>0 



1. (7) 



The last term on the r.h.s. of Eq (0) is a satisfactory discretization of the delta function provided the 
/fc's fall off sufficiently rapidly with k. The exclusion of /o ensures damping of transverse momentum 
integrals. In this approach, instead of a new induced quartic vertex we have introduced a short lived 
scalar, whose exchange simulates that vertex in a way that maintains Galilei invariance. Further, 
by leaving the choice of the /^'s open we might be able to tune their values to cancel unwanted 
symmetry violations induced by ultraviolet divergences in the continuum limit. The first term 
on the r.h.s. of Eq (^) is exactly what is needed to complete the modified cubic vertex when 
combining all of the contributions of the longitudinal gluons. 

Perhaps a more intuitive discretization would be to replace the derivative in D by a discrete 
difference: 

«-,Q.M»,-,.„, . (8) 

where the special treatment of the case k = 1 simulates the 6{x^) contribution we know must be 
there in the continuum. Unfortunately, this definition violates Galilei invariance because of the term 
that propagates only k — 1 steps in time: Newtonian mass conservation is temporarily violated. This 
causes considerable complications in calculations, but nonetheless displays interesting features. We 
will not pursue this option in the main text, but in an appendix we shall see that with a suitable 
counter-term Galilei invariance of the one loop self energy can be restored in the continuum limit. 



2.4 Tadpoles 

The exclusion of the propagators with zero r and zero p"*" renders every Feynman integral finite, 
making our discretization an effective regulator of divergences.^ But it also means that certain 
"tadpole" Feynman diagrams, which involve one or more propagators originating and terminating 
at the same vertex are excluded. In a theory with at most quartic vertices, these diagrams are 
limited to self-energy parts, which generally require counter-terms to enforce Lorentz invariance in 
the continuum limit. Thus errors induced by excluding these diagrams could simply be absorbed 
in the ultimate value of the counter-terms. 

Tadpoles arising from the cubic interaction (see Fig. |2|(a)) represent the vacuum expectation 
value of the color current density. The transverse components of these would vanish anyway because 
they are linear in transverse momentum but we can make sure all these tadpoles vanish by simply 
normal ordering the current density. However tadpoles arising from the quartic interaction, see 
Figj2|(b), would give a divergent non- vanishing result if the zero time propagator were inserted. As 
mentioned above one possibility is to absorb them in the self-energy counter-term. 

Another possibility is to note that, from the point of view of the continuum, one can just as well 
spread the quartic interaction over several time steps, in which case a candidate for the tadpole 

^The discrete light-cone quantization (DLCQ) industry which burgeoned in the mid eighties (for a review see Q) 
only exploited discretization leaving ultraviolet divergences unregulated. Discretization of ix'^ has the effect of 
introducing factors of e~°P (which is a popular way to regulate UV divergences) into loop integrals. 
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(a) (b) 

Figure 2: Tadpoles Feynman diagrams coming from cubic and quartic vertices. 




Figure 3: Quartic vertices get replaced by two cubic vertices and a fictitious scalar field. 

diagrams would emerge. We have already exploited this idea in our discretization of D , see Eq 
d?]) . A natural way to do this is to imagine that the quartic interaction is actually the concatenation 
of two cubic interactions mediated by a fictitious scalar field which is only allowed to propagate a 
few time steps^ We thus redraw the various quartic vertices as in Fig. ^. The fictitious scalars 
must be "ghosts" : to reproduce the quartic couplings, either their coupling to two transverse gluons 
must be taken imaginary or their propagator must be negative. We choose the second alternative 
for which the vertices are given in Fig. ^. Note that the quartic vertex which involves adjacent 




Figure 4: The couplings of the transverse components of the gauge field to the fictitous scalar field. 
The subscript m indicates that this is a magnetic ghost vertex. These replace the two quartic 
vertices at the bottom of Fig. |^. 

like direction spins in one channel can be viewed as a spin zero exchange in only one way, whereas 
the vertex with unlike adjacent spins in both channels becomes two exchange diagrams, giving a 
natural interpretation of the factor of 2 in the effective quartic vertex. If we now consider the 
"tadpoles" arising from connecting any pair of external lines, we see that there are 3 diagrams, see 

^Note that in higher dimensions, the fictitious field would be a transverse two-form instead of a scalar. Such an 
additional degree of freedom is presaged by the first order formulation of gauge theory in which Ff_ii, is treated as 
independent of and the Lagrangian density is -TrF'^/A + iTrf "(S^A^ - igA^,A^). Going to light -cone gauge m 

this formalism leaves, in addition to , the (nondynamical) fields and Fki. Our prescription simply gives these 

extra fields a short-lived dynamics. 
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Fig. ^, but the two with the topology of a cubic tadpole, which cannot be drawn in our discretized 
light-cone formalism, cancel. Thus the only remaining tadpole is the bubble diagram on the far 
right of Fig. H, which poses no problem for our formalism. 



o - -o - -o 



+ 



Figure 5: Three tadpole diagrams resulting from the spreading out of the quartic gauge vertex. 

The fictitious scalar propagator can be taken to be 

D{Q,M,k) = -/i.Toe-'^Q'/^MTo 



fc>0 



1, 



(9) 



where the /i^'s, like the /fc's, vanish rapidly with k. The normalization condition guarantees that 
the correct quartic vertex will be reproduced in the continuum limit. The exponential factor damps 
ultraviolet divergences, and with the form we have specified, maintains Galilei invariance even for 
finite lattice constants. Furthermore, by leaving the choice of the h^s open, we gain additional 
flexibility to cancel unwanted symmetry violations. The hope is that tuning the /^'s and /i^'s will 
remove the need for explicit counter-terms. 



3 Gluon Self-energy and Counter-terms at One Loop 

In this section we illustrate the way discretization regulates divergences by computing the one-loop 
contribution to the gluon self-energy part, iH^'^, defined as the sum of all one particle irreducible 
diagrams for the two-point function. For this purpose it is convenient to pass from representation 
to energy {E = p~) representation. With discretized r, this is accomplished by defining 

oo 

D'"'{Q,M,E) = ^e"'=-^Z)^^(Q,M,fe) (10) 

k=l 

oo 

U^"'{Q,M,E) = ^e"'=-^n^'^(Q,M,A:). (11) 



k=l 



The exact gluon propagator is then algebraically related to the the bare one and IT. We define 
g-Q /2MTo _ ^aE-Q /2MTo ^ terms of which the bare propagators have the values: 



u 



mQ.M.E) = (12) 



For D we have, for the discretization of Eq (0), 

I 

k>0 



Dr(Q-M,E) = |^-^-|L 53 (13) 
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Because of light-cone gauge, only IT'-' , H*"*" and 11"'"+ are required. By transverse rotational 
invariance, these quantities can be decomposed as 

ffi = Q^Qmi + r]'m2 

ff+ = MQm[ 

U++ = M^U'I (14) 

We shall find that at the one-loop approximation, 11'/ = 11'^ = Hi and if that were to hold generally, 
the exact propagators would be given by 



2M1-U- UU2/2M 
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2M2 l-u- UU2/2M 



b- = ^ 



3.1 Calculation of Hi 



+ 



1 + ToBi X^fc /fcii'' 2MTol-u-uYl2/2M 



(15) 




Figure 6: The two bubble diagrams that contribute to 11"'"+. 

After these preliminaries, let us now turn to the computation of IIi at one loop. The simplest 
term is H"'""^, which is given in representation by (note that there are 2 equal diagrams that 
contribute (see Fig. P)) 



-e 



-fe(pV2i+(Q-p)V2(M-0)/To 



_ 2g^N, 1 feQV2MTo N^aHM - 1)(M - 2) ,Qy2Mn (.^^ 

^ 1=1 ^ 
= ^|^(M-l)(M-2)(-ln(l-^)). (17) 



Figure 7: The four bubble diagrams that contribute to 11 
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The evaluation of II-'"'" is not much harder. Here there are four diagrams that differ only in 
the prefactor of transverse momentum (see Fig. 0). Since the prefactor is linear in momentum, we 
need only to remember that the Gaussian integral over transverse momentum involves the shift 
p — > p + IQl/M. After this shift the term linear in p integrates to zero, so the net effect is to set 
p = ZQ/M. Thus the net prefactor is 



^{21 - M)[(l + l/M) + (-2 + l/M) + r,+_(-l - l/M) + r?+_(2 - l/M)] = 



which involves the identical sum over / as 11^+. Thus we end up with 











Figure 8: The eight bubble diagrams that contribute to 11^'^. 

The coefficient of Q^Q^ in II*-' may be singled out by computing 11^^. This leads to a similar 
calculation to the above because, after the shift in integration variable, and p^p^ both integrate 
to zero. Thus again, prefactors of p^ may simply be set to IQ^ /M. Eight diagrams contribute to 
this quantity, equal in pairs (see Fig. |8|). The prefactors of these diagrams combine as follows: 

^ y ((M + If + (2M - If - (M + /)(2M - /) - (M + 1){2M - I)) = ^ ? {21 - Mf. (20) 
So again we have the same sum over /, leading to: 

The upshot of the calculations so far is that 

The equality of the various H's holds even at finite m, a. This can be understood because the 
diagrams we have evaluated show no violation of Galilei invariance. Our result for Hi can be 
compared with the result from the study of asymptotic freedom in the infinite momentum frame 
i: 
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in which a simple cutoff < was employed. To make the comparison, note that in the 
continuum limit, a ^ 0,M ^ oo, 1 - u (Q^ - 2Q+E)/2MTo = Q"^ /2MTq. Thus for us, the 
role of A-^ is played by the quantity 2MTq = 2Q^ /a. Note that if we choose to keep Tq fixed, the 
ultraviolet cutoff is removed by simply taking M ^ oo. 



3.2 Calculation of Ha 

Finally we turn to the real core of the self energy calculation, the determination of 112, which can 
be inferred from = Q^Q'^Iii + n2 = n2 + Q'^Iii/2. Seventeen dia grams contribute to this 
quantity, fifteen of which do not involve D and are relatively simple to analyze. First consider 




Figure 9: The three bubble diagrams that contribute to 11^ involving only transverse components. 

the three diagrams that only involve transverse internal propagators (see Fig. P). The prefactors 
of transverse momentum combine as follows: 



{2p - Q)\2p - QY + {p + Q)\p + QY + (2Q - pY{2Q - pf 

^ 3(p2 + (1 - l/M + (//M)2)q2), 



3(p2 + Q2 _ p . Q) 



(24) 



where in the last line we have indicated the result of shifting p by IQ/M and dropping the term 
linear in p. The transverse integral of the term in is the same as before, but the term in p^ 
gives an extra factor 2TqI{M — l)/kM. These terms require the sums 



M-l 



(25) 
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1=1 



respectively. The contribution of these three diagrams, , is given by 
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The next class of diagrams consists of the eight graphs with one D^^ propagator, shown in 
Fig. Remembering the single factor r]-^ , we find that the prefactors combine as follows: 
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Figure 10: The eight bubble diagrams that contribute to Hg . 



where as before the arrow indicates the effect on the prefactors after the usual shift of integration 
variables. Performing the by now familiar integrals and sums leads to the following result for the 

tAV 



contribution of this class of eight diagrams, labeled by Ut'^- 
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Combining this with the previous three diagrams yields for the eleven diagrams considered thus far 
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Figure 11: The two bubble diagrams that contribute to Hq^. 

The prefactors in the pair of diagrams (shown in Fig. [Tl| ) involving two D^^ propagators combine 
to form 



(M + /)(/-2M) 



p(Q-p) 



(M + /)(2M-/) 2 
1{M - I) ^ 



(M + /)(2M-0 
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This time the required sum is just 

M-l 

{2M^ + IM - f) = M{M - 1)(13M + l)/6, 



M2 



(30) 



(31) 



1=1 



so that these two diagrams simply double the result of the first eleven. So, in summary, the total 
contributions of the thirteen diagrams that do not involve D to 11^^ and IT^^ are 

g'^Nc 13M2 - 12M - 1 [ 1 " 
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(32) 
(33) 
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Figure 12: The two bubble diagrams that contribute to which involve a D propagator. 



The two diagrams in Fig. 12, have a D propagator as one of the internal lines and lead to a 
qualitatively different evaluation than the first thirteen. Using (^) the Gaussian exponents are the 
same as with the first thirteen diagrams. One finds 
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where we have made use of the identity 
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(35) 



where ^{z) = V {z) /T{z) is the digamma function and 7 is Euler's constant. At large M, we have 

V.(M)~lnM--i--5],^, (36) 



2M ^ 2nM2'^ ' 

n>l 



where B2n are the Bernoulli numbers. 



Figure 13: The two bubble diagrams that contribute to which correspond to the tadpole 
contributions induced by the spread out quartic vertices. 

With the quartic vertices realized as the exchange of a magnetic scalar, we have a definite 



each give equal contributions. Calling the tadpole contribution , we find 



proposal for the tadpole contribution to the self-energy, namely the two diagrams of Fig. |13|, which 

tAV 

O.r M-l 



'(27r)3T( 



hk 



-k[Mp^ /l{M~-l)+Q^ /M]/2To 



2{M -I) 



2 ,r M-l 

g \ - '"fc A:Q2/2MTo 



E 



4vr2 ^ Mk' 
1=1 

87r2 ^'"^ k 



(37) 
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Adding + 11^^ to the first thirteen diagrams, gives (in energy representation) 

nAv ^ ftAv + ^ jjAv ^ nAv ^ nAv 



47r2 



\2t 



(38) 



To get II2 we must subtract Q ni/2, which exactly cancels the last term leaving: 



-.2 AT r °° „,fc 



g|!(.Mi.(M).,- <^-^Kn^-^' ) 

A:=l A;=l ^ ^ 



(39) 



The first term inside the large square brackets is exactly the contribution one would obtain using 
the modified cubic vertex (^) in the three diagrams of Fig. |9|. The terms with represent the 
tadpole diagram that would have been constructed from the induced quartic vertex had been 
eliminated from the formalism. This is as it should be because summing over all the diagrams 
involving and D is the graphical equivalent of the elimination of the j4+ degree of freedom. 
The formalism with eliminated gives an extremely efficient calculation of the three bubble 
diagrams, as shown in the next subsection, but no candidate for the tadpoles. But by treating Aj^ 
graphically we are led to a useful proposal for the tadpole diagram. 

In order to compare our results to the continuum calculation of Ref. we must examine the 
limits It — > 1 and M ^ 00. The behavior of the first term of Eq. (|39| ) in this limit is transparent 
once we use the identity 

Y^^ = ^ \nu In(l-n)- j;^^^^^ + (l-^)(ln(l-n)-l), (40) 
fe=l k=l 

for « « 1. The last two terms can be expanded about u = 1: 

k=l k=l k=l 

00 , 

= E^ + {u-l) + 0{{n-in (41) 

k=l 

and a similar expansion for f^u^ /k. The coefficient of the linear term is completely determined 
by the normalization conditions on the /'s and /I's. Momentum independent terms in 1X2 imply a 
(divergent and noncovariant) tachyonic gluon mass squared in perturbation theory. Since the gluon 
has only helicity ±1, such a mass is inconsistent with Poincare invariance of the continuum limit. 
Clearly these symmetry violating terms would be cancelled if we could impose the constraints 



00 „ 2 



k=l 



k=\ 
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but clearly the second of these is impossible with hk independent of M (implicitly assumed in 
d)). The best we can do is to set the r.h.s. of Eq. (|4^ to — vr^/lS, which would cancel the linear 
divergence of 112. However, there remains a finite gluon mass which is still inconsistent with the 
Poincare invariance of the continuum theory. One option would be to cancel this with a mass 
counter-term, whose coefficient would have to be determined order by order. 

Another approach is to allow hk to be dependent on M. We have the freedom to make this 
replacement as long as we recover the correct continuum limit (M oo and a — > 0) of the theory]^ 
With this modification, Eq. (|37|) is replaced by 
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where the sum over I is not performed. Similarly Eq.'s (| 
Thus 112 is given by 
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(44) 



should now include the 



corrected 11^^ 
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and the expansion of Eq. 
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(45) 
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The constraint equation that replaces (|4^) is 
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(46) 
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which can be satisfied if we require 
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k=l 



hk{l) 
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1 



(47) 



(48) 



We prefer this approach to that of a mass counter-term, since it is possible that this is a uniform 
description that works order by order (at each order the cancellation places constraints on higher 
moments of h). The hope is that this can be used non-perturbatively. 

In the continuum limit M(l — u) — > Q'^ /2Tq stays finite and we find that 112 tends to 



Ho 



167r2 To 



8(lnM + 7) 



22 
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In 



2MTq 3 



(49) 



Remembering that our 112 is a factor of —Q'^/Tq times that defined in Ref. |p, we find agreement 
for the coefficient of InQ^^ provided we identify Me^ with Q'^ /e and 2MTq with A^. We do not 
get, nor should we expect, the same (finite) coefficient of Q^. 

^We have the same flexibiUty for the /'s in (j^ but as we will see this is not necessary at least at one loop order. 
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3.3 A Brief Calculation of U 



In the work just completed, we deliberately kept the degree of freedom in the graphical rules 
in order to keep the calculation as close as possible to one in other gauges. However, having seen 
how all of the graphs with longitudinal gluons combine so nicely, it is appropriate to note that the 
calculation with explicitly eliminated, so that the Feynman rules refer only to the transverse 
gluons, is much more compact and efficient. With our prescription we can exploit this simplification 
by using the modified cubic vertices Eq. (^) for the transverse gluons. At the same time we retain 
our replacement of both the bare and induced quartic vertices by the exchange of two short-lived 
scalars. The self-energy diagrams involving those scalars will be exactly as described previously (i.e. 
the terms involving /'s and /I's). However, all of the remaining contributions to H^' are reduced to 
the two or three diagrams involving only transverse gluons and the modified cubic vertices. 

Only two diagrams contribute to H^^, the last two diagrams in Fig. ^, and they each involve a 
prefactor 



{Mp - lQ)^{Mp - IQ)^ M'^p^p^ 



(50) 



after the shift in momentum. Clearly this integrates to zero so H^^ = 0. Finally there are only 
three non-tadpole graphs to consider for H^"^, namely those in Fig. In this case the relevant 
prefactors from the three diagrams contributing to H^^ combine as 



P{M-lf 
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(M - /)2M2 
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(M - ly 

- M2 



+ 



(Mp - IQY 



+ 



(M - If 



P{M-lf {M-lfM^ M2;2 



M2p2. (51) 



after the usual shift in momentum. Note that because these vertices are manifestly Galilei invariant, 
there is no term proportional to in the prefactor. Thus after integration over loop momentum 
we are left with the contribution to H^"^: 
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(52) 



In this way of organizing the calculation, the longitudinal components of H'^'^ play no role and the 
new diagrams contributing to H^^ give Yi2 directly. For the non-tadpole part we find: 



H; 



non~tadpolc 



47r2 



E> 



k=l 



(53) 



which is clearly a much simpler and more compact calculation! To this result must be added the 
contribution of the fictitious scalar diagrams, which represent the tadpoles. Needless to say, the 
calculation of more complicated processes should make use of these new Feynman rules, which we 
have summarized in Fig. 14. 
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Figure 14: Summary of discretized Feynman rules using only cubic vertices. We have explicitly 
inserted a factor of 1 /Tq for each vertex arising from the discretization. 

4 Summing Planar Diagrams: Fishnets 

The discretized Feynman rules given at the end of Sec. ^ provide a tool to sum classes of diagrams. 
As described in Q summing diagrams on a light front has a direct interpretation as the path 
history quantum evolution of a system of particles moving in the transverse space under Newtonian 
dynamics. By fixing the total discretized = Mm, the maximum number of particles present 
at any time is M. Because the vertices allow particles to fuse and fission, particle number is not 
conserved and there is quantum mechanical mixing between states with any number of particles 
between 1 and M. 

We are particularly interested in the class of planar diagrams singled out by 't Hooft's Nc ^ oo 
limit of QCD. It is actually more precise to think of this class of diagrams as drawn on a cylinder 
rather than a plane: At any time the system of particles is ordered around a ring, and interactions 
only exist between neighbors on this ring. Thus the stage is set for the particles to bind into a 
closed polymer chain. This was previously investigated in where scalar matrix field theory with 
quartic couplings, ATr(/)^/4 was considered. These interactions are attractive (repulsive) if A < 
(A > 0). Thus bound chains can form only if A < 0, the unstable sign. In the interpretation of 
the sum of diagrams as a sum over histories of a system of particles, this sign assures that all 
histories contribute with a positive weight. Using the discretization of and ix~^, as reviewed 
here in Sec. |2|, a strong 't Hooft coupling limit AA'^c oo was formulated and analyzed. This 
was achieved by focusing attention on the cylinder diagrams that evolve a system of particles with 
P^ = Mm a fixed large number N steps forward in time. The limit singles out those diagrams in 
which every particle has the minimum P^ = m (so the number of particles is maximal = M), and 
each propagator evolves only one time step. For even M these diagrams include the large fishnet 



16 



diagrams that form a seamless web of quartic vertices and propagators, and the resulting Gaussian 
integrals could be identified as a discretized path integral for a relativistic bosonic string quantized 
on the light-cone. 

This work on scalar field theory was immediately followed by a first attempt to apply these 
ideas to large Nc QCD In that work QCD was formulated on a light front, with P+ and 
ix~^ discretized. The ordinary bare vertices of QCD, both cubic and quartic, were used. Because 
the quartic coupling in QCD is of order g^, the literal strong coupling limit, as formulated in Q, 
favored fishnet diagrams with only the primitive (i.e. non-induced) quartic couplings. Actually 
this conclusion required the ad hoc exclusion of the = exchange part of the induced quartic 
interactions arising from fixing the light-cone gauge. Nonetheless, the resulting fishnet was very 
interesting: the spin of the gluons played the role of the arrows of a certain six vertex model, known 
as the F- model In fact the four gluon vertices of the field theory were exactly the vertices of 
the F-model. The fact that some of the vertex weights were negative did not cause problems for 
the leading strong coupling fishnets because those diagrams always had an even number of negative 
weight vertices. However the problem with them reappears at next order because the deletion of a 
single negative vertex reveals a repulsive nearest neighbor interaction in that spin channel^. Thus 
that channel could not have formed a bond in the first place. In some spin channels there were 
also positive weights, so that the bonds could form. However, unfortunately for these fishnets, 
the attractive channels are ferromagnetic: the only long polymers that could be formed by these 
interactions would have enormous spin! 

The problem is that at strong coupling only the quartic interactions survived with the dis- 
cretization of The spin-spin interaction from gluon exchange has anti- ferromagnetic behavior, 
and it is possible that a discretization that allowed the exchange interaction to compete with the 
quartic interaction could cure this problem. To explore this possibility, one of us examined the 
relative strengths of quartic and cubic exchange interactions for neighbors on a gluonic chain by 
putting the two gluons in a spherical MIT bag [10|. In that context one can see explicitly, not only 
that the cubic exchange of a transverse gluon is anti-ferromagnetic, but that its strength (at least in 
weak coupling perturbation theory) is more than sufficient to reverse the ferromagnetic character of 
the quartic interaction. A major shortcoming of the discretization of [|| is that at strong coupling 
the cubic interactions have no opportunity to compete with the quartic interactions. 

The discretization developed in Sec. § of this paper is more promising. In fact, we do away 
with quartic interactions completely: All interactions are cubic! The quartic interactions have 
been replaced by the exchange of fictitious scalars, and these exchanges are not enhanced by strong 
coupling over the exchanges of the ordinary transverse gluons. Without quartic vertices the basic 
cells of the densest diagrams are no longer square but hexagonal: the fishnet looks like a honeycomb 
(see Fig. 15). These fishnet diagrams require particles with both 1 and 2 units of P^. Thus the 
leading fishnet structure is somewhat looser than in the quartic coupling case. The inclusion of 
two different values of P"*" in the leading approximation also sets the stage for the emergence of 
a string degree of freedom (provided string states do form!) corresponding to fluctuations in P^ . 
Such a degree of freedom is expected to be described at long worldsheet wavelength by Polyakov's 
Liouville field lll|. 

Because all of our vertices are cubic, the paradigm scalar field theory is now gTj:(p^/3 (see 111 
for a discussion of this model in 1 + 1 dimensions). Of course, the presence of factors of transverse 
momentum in the QCD cubic vertices will cause a profound qualitative difference between gauge 
and scalar theory. For one thing, the scalar theory is super-renormalizable with the cubic coupling 

^ Contrast this with the cubic scalar theory where the sign of the coupling is indeed irrelevant because the deletion 
of a single cubic coupling is not allowed by the Feynman rules. 
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Figure 15: Hexagonal cell fishnet made of only cubic vertices. The dashed lines indicate the closed 
cylindrical topology. 

carrying dimensions of mass. This means that weak or strong coupling is determined by the size 
of the ratio g^Nc/fi^, where /i is a mass scale relevant to the calculated physical quantity. For 
example, in our strong coupling considerations, we can simply take = Tq. Also, since (j)^ is 
unbounded at large (p in one direction, the theory is ultimately unstable, although this instability 
is not evident in weak coupling perturbation theory if the scalar field has a non-zero mass. But the 
topology of graphs and momentum flow is the same in both scalar field theory and our version of 
discretized QCD, so we shall exploit the scalar paradigm to illustrate such common features. The 
absence of spin degrees of freedom in the scalar theory is a helpful simplification for at least some 
issues. 

Before turning to details, we comment on the ambiguities in our setup contained in the values of 
the /'s and /I's (see Eqs (j^ and (^)). As shown in our study of the gluon self-energy in Sec. ^, weak 
coupling perturbation theory constrains moments of these quantities (for example, at one loop we 
find constraints on ^ fk/k and Yl^k/k). Since strong coupling emphasizes short times, we expect 
this limit to put constraints on the fk and for low values of k. Thus we gain complementary 
or "dual" information about the theory as we explore both limits. It should be stressed here that 
there is no "unique" discretized theory to associate with continuum QCD: all sorts of lattice details 
get washed out in the continuum limit. Instead our goal is to find a single lattice model that shows 
continuum QCD at weak coupling and string theory at strong coupling. This dual requirement will 
hopefully help to determine a unique theory. 

As in iQ] we shall consider the sum of cylinder diagrams that evolve a system with P+ = Mm 
forward in time by the amount T = Na. For fixed initial and final states, according to our 
prescription, there are only a finite number of diagrams that contribute. This defines a definite 
model that can be studied as a function of the bare couplings. If the complete sum could be done for 
arbitrary M and A^, one could then read off the exact spectrum of the continuum theory by studying 
the limit M,N oo, with parameters tuned so that the limit is nontrivial. This is presumably 
too ambitious, but in the next section we shall at least be able to deal nonperturbatively with 
some small values of M. One might also envision studying moderate values of M numerically on a 
computer. In the rest of this section we shall discuss the fishnet diagrams that describe the infinite 
coupling limit of our model. 

Let us first consider the scalar paradigm. It is sometimes helpful to define a transfer matrix T 
which evolves one step forward in time. In order to do this, start with time continuous and express 
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the exact time evolution by an amount a in the interaction picture: 

This expression is exact, and of course it does not correspond to any discretization. Our dis- 
cretization is given by approximating each Vi{ti) ~ ^(0) = V and only retaining the term 
T[Vi{ti) ■ ■ -Vi^tn)] when each Vj acts on a different subsystem of the particles present initially. 
We shall therefore write the transfer matrix for our discretized theory as 



oo 



r = e-"^« ^ ^-j- [y"] , (55) 



n 

n=0 



where we understand [V^] = unless each of the nV^s acts on a different subsystem of the particles 
present. With this understanding we implement our discretization rule that every line in a diagram 
propagate at least one step in time. These approximations are strictly valid for sufficiently small a 
at fixed coupling parameters. But we use Eq. ^ to define a discretized fishnet model at fixed finite 
a, which we intend to study at all values of the coupling, including g — > oo. Although the strong 
coupling limit at fixed a (as always) takes one far from the original continuum theory quantitatively, 
we hope that it will lead to a new continuum QCD string theory bearing qualitative resemblance 
to real QCD. But there is, of course, no a priori guarantee of this outcome. 

An efficient way to implement the Nc ^ oo limit is via the Fock space approach of One 
chooses a state of the form 



M 
i=l 



\j^Y. /Tr[ai,^(pi)---aj^^(p^)]|0)^,(Pi,^i,---,P^,^^), (56) 



and applies the transfer matrix keeping only terms that survive the Nc ^ co limit. The second 
sum is over all partitions of M, such that Mi + M2 + • • • + -M^ = M. All such terms retain the 
color trace structure and describe interactions between nearest neighbors as defined by the color 
trace. If we wish also to take the infinite 't Hooft coupling described by the densest fishnet, we 
choose M even (M = 2n) and restrict £ in the sum to the two values: i = 2n with = 1 for all 
k and i = n with Mk = 2 for all k, and we require that every particle present participate in an 
interaction at each time step. This leads to the following coupled equations for and ip2n when 
is an eigenstate of the transfer matrix: 

gnAr^"/2g-E,(P?+M?,)/2To ^ 
tV'2n(Pl,---,P2n) = (4To)"(27r)3"/2 2 "^P^ ^ , PS + P4, " " " , P2n-1 + P2n) 

+ V'n(P2 + P3, P4 + P5, • • • , P2n + Pi)] (57) 

2<7"Are"/'e-S^('i?+^o)/^^" frr 2 
tVn(qi,---,qn) = (4To)"(27r)3"/2 y J_J_d A;j'i/'2n(ki,qi -ki,---,kn,qn-k„),(58) 

where t is the eigenvalue of the transfer matrix, and we have included a bare mass fiQ for the scalar 
field. In the continuum limit t = e~"^. We have suppressed the M^'s in the arguments of ipn and 
tp2n due to their simplicity. Clearly, we can eliminate ■02n to obtain a single equation for ^pn■ 

tVn(qi,---,qn) = (^ \6ro(2vr)3 J « ^ 11 ^[^'^(qi, q2, • • • , qn) 

+ M^i + k2 - ki, • • • , q„ + ki - k„)]e-S.[^?+('i'-l'>)']/2'^°. (59) 



19 



We see that the /xq dependence is a trivial factor in this strong couphng equation, so we shah set 
yUo = in the following. The integral of the first term in square brackets is elementary yielding a 
factor Tr"'e~^i'^i^^'^° . Defining 

and rearranging the equation leads to 

j 

= X-e~^.^V^^^ J n^V'n(q./2 + qm/2 + k.+i-k,)e-S^'^?/^o, (60) 
j 

where in the last line we have shifted integration variables to complete the square in the Gaussian 
exponent. This equation sums diagrams including not only the basic fishnet, but also fishnet 
diagrams containing any number of time intervals in which n subsystems each with M = 2 propagate 
freely for arbitrary lengths of time (see Fig. |l^). This complication is described by the second term 



Figure 16: Fishnet with a section of M = 2 subsystems propagating freely. The dashed lines again 
indicate the closed cylindrical topology. 

on the l.h.s. In considering the effect of this term, one should keep in mind that no self mass counter- 
terms have been included in the derivation of (|60|). For example, in weak coupling perturbation 
theory, the self energy bubble by itself would lower the scalar mass squared by an infinite amount. 
To keep the scalar mass non-tachyonic at weak coupling one needs a mass counter-term. With 
discrete time, it is convenient (as in our treatment of tadpoles) to spread such a mass counter-term 
over several time steps by introducing a short lived fictitious scalar with a quadratic coupling to 
the real scalar field of order g. We would then have additional diagrams as in Fig. and the 



Figure 17: Additional diagrams introduced by spreading out a mass counter-term, 
upshot for Eq. (^) would be an adjustable coefficient in front of the second term on the l.h.s.: 
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/ n + ^^+1/2 + ^^+1 - k^)e-^«''»'/^^ (61) 



We shall begin by dropping this term completely {i.e. by choosing 61 = 1) to determine the 
contribution of the basic fishnet for the scalar cubic theory. Later, we shall comment on the effect 
of 61 < 1. 

In order to understand the dynamics inherent in the basic fishnet, it is helpful to express the 
integral transform on the r.h.s. of ( |60|) as an operator in the state space of a first quantized system 
of n particles. It is straightforward to show that the appropriate operator is given by 



(62) 



where the p's and x's are the momentum and coordinate operators for the n particle system. Here 
$7 is an operator defined in momentum and coordinate bases by 



^ = / Ild^Pj\Pk){MkiPi\ = j \{(fxj\iLiMik){-i^k\ 



and M is an n X n matrix defined by 



MkiPi 



Pfc + Pfc+i 



(63) 



(64) 



It is easy to check that M is invertible provided n is odd. For n even there is a zero eigenvalue 
which must be separated and handled before continuing the analysis. For simplicity, we assume n 
is odd in the following discussion. One can readily verify that $7 has the following action on the 
coordinates and momenta: 



Mkipi^l. 



(65) 
(66) 



Because of the Gaussian structure of T, it also has a linear action on the coordinates and 
momenta: 



Txfc = [xj+ipz/ro](M-i)ifcT 

TP. ^ (m.p, - . *,/T„,[2(M-.)„ - (M-.)„« - (M-.)„._j) T. 

This linear action can be diagonalized by passing to normal modes: 
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One then finds that the modes I all decouple from one another under the action of T: 
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(72) 
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We can now search for eigenoperators of the form + ^/p/. This leads to a quadratic equation 
for £,1, 



1 



n] 



0, 



with solutions 



^'^"4ro (^^\/^+sin2(/vr/n) 
These eigenoperators change the eigenvalue of T by a factor 
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(74) 



Note that these eigenvalues are not real because T is not a Hermitian operator. However, also 
note that the second factor is positive for both branches and for < / < n. The eigenvalue is 
therefore always in the right half complex plane. We also have that |A_(_A_| = 1 which implies that 
|A+| > 1 > |A_|. Moreover, the first factor which contains the complex phase can be rewritten in 
two ways 
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(75) 



which shows that the phase is proportional to the fishnet momentum created by the eigenoperator: 
l/n for I < n/2 and — (n — l)/n for I > n/2. Cyclic symmetry of the initial wavefunction implies 
that the total fishnet momentum must be 0. 

The ground state (belonging to the largest eigenvalue T) is determined by the condition that it 
be annihilated by all the eigenoperators which increase i^, x; + £,^pi- Its wavefunction is therefore 
proportional to the Gaussian (with normalization M) 
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(76) 
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which is always damped because x^* = x„_/. The eigenvalue corresponding to this state can be 
obtained in the following way. First, we observe that 
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Then together with Eq. (|6^) we get that 
T^G = AA22("-i)exp 
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Thus with the use of the identities 
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and 
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we find that the eigenvalue for the ground state is given by 

(n~l)/2 



n 

1=1 



1 



1 



IH — sin (/vr/n) H — sin(Z7r/n)A/8 + sin {lir/n 



(81) 



Since it is positive, all cyclically symmetric states generated by applying suitable monomials of the 
eigenoperators to the ground state will have positive eigenvalues of T. 

Clearly the long fishnet wavelength excitations show behavior identical to those of the continuous 
light-cone quantized bosonic string. The excited states are obtained by applying appropriate zero 
momentum monomials of the eigenoperators + S,fpi to ^g- From the interpretation t - 
we see that 
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where we have used the Euler-Maclaurin summation formula for large n 
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In lattice string theory the bulk term proportional to n contains no physics and can be dropped 
(see |]l4l). The ground state string mass squared is predicted to be (recall that P"^ = 2nm) 
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We also see that the basic energy splittings are given for / ^ n by 
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or splittings in mass squared of 
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This result shows that the string arising from our basic cubic fishnet has an effective rest tension of 
To/\/2, corresponding to a Regge slope parameter a' = I/^ttTq. Noting that here the transverse 
dimensionality is 2 and not 24, Eq. (84) gives the usual result of bosonic string theory, —d/Qa'. 

This is all for the basic cubic scalar fishnet, in which the second term on the l.h.s. of Eq. ( |60|) 
is tuned to zero. Including that term, we find a solution for general n intractable. However, quali- 
tatively, we can say that it introduces a continuum threshold at i = A"'/^(l — (Ji)"/^, corresponding 
to a threshold energy Eth = — (n/2a) ln(A(l — 5i)). As long as the basic fishnet described above 
produces a ground state energy Eg < Eth-, we can expect qualitatively similar physics for large n. 
However, for Eg = Eth we can expect that the seamless fishnet structure begins to be disrupted 
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Figure 18: Typical example of a strong coupling QCD fishnet. 



with a dramatic qualitative change in the physics. We shall explore this effect for small values of 
M in the next section. 

Finally, let us turn to QCD. The first major difference is that each line can exist in four different 
internal states, corresponding to the two polarizations of the transverse gluon (with spin = ±1) 
and the two fictitious scalars introduced to simulate the quartic interactions (see Fig. 15 for an 
example of a generic dense QCD fishnet lattice). In Fock space language this can be described by 
affixing a 4-valued index to the creation operators. The basic fishnet diagrams will be the same as 
in the cubic scalar theory, but with the complication that the vertex value depends on the states 
of the lines entering it. In particular, some of the vertices are linear in the transverse momenta of 
the incoming lines, leading to an interesting spin-orbit coupling on the fishnet world sheet. Thus 
the QCD fishnet dynamics requires the solution of a two-dimensional lattice spin system with a 
nontrivial interplay between the internal "spin" variables and the structure of the lattice itself. 
Contrast this with the fishnet contemplated in and based solely on the quartic coupling. In the 
latter situation the spin degrees of freedom decoupled from the transverse coordinate degrees of 
freedom and corresponded to the soluble F-model, one of the 6-vertex models. The fishnet model 
we are proposing here has a considerably richer structure, which we shall begin to explore for small 
values of M in the next section. 



5 Baby Fishnets 

In Section |^ we discussed how our discretized Feynman rules can be used to determine the dense 
QCD fishnets for strong coupling. This was also discussed in the context of the paradigm cubic 
scalar theory. In this section we are interested in the dynamics of the discretized theory away 
from strong coupling. Potentially the strong coupling limit will not be described only by the dense 
fishnet lattice. While the ultimate goal is to do this for M ^ oo (remember = Mm), we will 
begin by analyzing systems with small values of M (i.e. baby fishnets). 

The simplest non-trivial QCD fishnet has M = 2. We always understand our fishnets to 
propagate color singlet systems so that they have cylindrical topology. Then the M = 2 fishnet 
has no interesting dynamics due to the fact that color-singlet gluons decouple from gluon bubbles. 
However, the color adjoint M = 2 gluon propagator which plays an important role as a subsystem 
of larger diagrams, can be solved to all orders in perturbation theory due its simplicity. We leave 
investigation of M > 3 to the future. Another possible avenue of investigation is systems involving 
sources rather than pure glue, (such configurations are discussed in |l^). We also defer exploration 
of such systems. 
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In Section ^ we were only able to solve the strong coupling cubic scalar fishnet for general n 
with the term on the l.h.s. of Eq. (61) cancelled via a mass counter-term (we will define = l — 6i). 



But for special cases of M = 2n we can solve (61) including K^-term. The simplest of these is the 
M = 2 scalar fishnet, however, in this case the only effect of the «;^-term is to rescale the solution 
presented in Section ^ For the more interesting cases of M = 4, 6 we will see that they too can be 
solved. For M = 6 we restrict attention to the s-wave sector. 

5.1 M = 2 States of QCD 

The M = 2 color singlet glueball states display no dynamics, because our discretization with exclu- 
sively cubic vertices only allows interactions via mixing between one gluon and two gluon states, 
and there is no interacting color singlet gluon. (Even if the gauge group is U{Nc) the abelian gluon 
completely decouples in the pure gauge theory.) To understand this decoupling in terms of our 
Feynman rules, note that on a cylinder the gluon self energy bubble can close in two ways as in 



Fig. |19[ It is then easy to see from our rules that the two diagrams are equal in magnitude and 
opposite in sign. The conclusion is that the M = 2 color singlet channel consists of two free M = l 
gluons or, in the U{Nc) case, a single free M = 2 gluon. This trivial situation is due to the manner 
in which the quartic vertices of the initial gauge theory have been replaced by scalar exchange. 
Only the "direct channel" scalar exchange is allowed at M = 2, and the part of the quartic vertex 
that is described by the "crossed" channel exchange only makes its appearance for M > 3. 








(a) (b) 

Figure 19: Here we see how bubbles of a generic M = 2 fishnet lattice can close in 2 ways due to 
the cylindrical topology. The wraparound diagram (a) has the opposite sign to the bubble diagram 
(b), so that they cancel. 

Thus the only nontrivial M = 2 channels are color non-singlets. Moreover, fixing M = 2 limits 
the allowed diagrams so drastically that the nonvanishing ones can be explicitly summed to all 
orders in perturbation theory. We first look at the M = 2 gluon propagator, which can be simply 
read off from Sec ^. For simplicity we work in the center of mass frame. The diagrams that 
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Figure 20: Gluon propagator for M = 2. 

contribute to the gluon propagator are depicted in Fig. |2^. The shaded bubble corresponds to all 
the one-loop bubble diagrams that contribute to the transverse gluon self energy, 11^^ which is 
obtained by putting Af = 2 in Eq. 
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where u 



in the center of mass and with 
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These constraints on the /'s and /i(l)'s have been determined in Section ^ at large M (see ( ^2|) 
and (^)) in order to cancel divergences in 112. We tentatively impose the same constraints at all 
finite M in order to have a uniform description for all M. The exact transverse gluon propagator 
for M = 2 is (see ([|)) 
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(89) 



This propagator evolves a spin 1 color adjoint system, which by itself would not correspond to a 
glueball, which must be a color singlet. Because of its importance for larger diagrams, it is worth 
understanding the energy eigenstates implied by the propagator's pole structure. The (1 — u)~^ 
factor out front is just the massless gluon pole {E = implies u = 1). 

Zeroes of the quantity in square brackets in Eq pi determine any additional eigenvalues: 
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(90) 



Even with the /'s and /i(l)'s general, one can note that the r.h.s. tends to — oo as u — > 1, with 
behavior completely fixed by the constraints. Also for u — > the r.h.s. vanishes quadratically as 
(9 — 9/i — /ii(l))u2/4. Therefore if 9 — 9/i — hi{l) > there would be at least two solutions 
for sufficiently large g, the lowest of which would tend to n = as (7 — > 00. If, however, the 
inequality were reversed, the r.h.s. might never be positive in which case there would be no solution. 
Alternatively, if it did cross the axis there would be at least two solutions the lowest of which would 
tend to some nonzero n = > as 5 — > 00. It is amusing to see which of these behaviors is 
suggested by a minimal solution of the constraints so far imposed (see (^8|)). For M = 2 we can 
meet the constraints with only the first two elements of each series nonzero, which are then fully 
determined by the constraints: 
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The eigenvalue equation then reduces to 
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(92) 



where Li2(ii) is the dilogarithm (Spence) function [16|. As we can see in Fig. ^ this minimal 
choice shows no physical eigenvalue, since there is no positive solution for 8tt'^ /g'^Nc for any value 
ofue [0,1]. 

Finally, the M = 2 color adjoint magnetic scalar propagator also receives self energy corrections 
which can be summed exactly (see Fig. [2^ ). (Note that the fictitious electric scalar (solid line 
propagator) does not play a role in the M = 2 channel, since it's coupling to two M = l transverse 
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Figure 21: A plot of Eq. ( |9^ ) with S-n"^ / g^Ng. (along the vertical axis) against u. 




Figure 22: Magnetic scalar propagator for M = 2. In each bubble the transverse gluon index may 
circulate in either direction. 

gluons is zero.) Although the magnetic scalar's contribution to the dynamics of a color singlet 
M = 2 glueball cancels, its propagator describes a spin Af = 2 color adjoint subsystem in larger 
fishnets, and so it is also useful to analyze it here. In these diagrams the bubbles correspond to the 
one-loop self energy diagrams of the fictitious magnetic scalar (the dashed scalar). The magnetic 
scalar self energy (see Fig. p3| ) is given by 
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The bare magnetic scalar propagator for M = 2 is then 
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where the /ifc(2) have to obey the constraints 
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The exact propagator is then given by the geometric series 




(96) 



27 



Figure 23: The two bubble diagrams that contribute to the magnetic scalar self energy, n^- 

We again investigate possible energy eigenstates by looking at the pole structure of this amplitude. 
Focusing on the denominator we see 
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Again we see the same possible behaviors as in the case of the gluon propagator (except, of course, 
there is no massless pole at u = l!). In this case we also present the results of choosing a minimal 
set of the /ifc(2) to satisfy the constraints. Doing this yields 
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With this set of parameters the denominator factor in (97) will have a pole if there is a solution to 
the following equation 
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where a = /i2(2)//ii(2) w —1.646. As we can see in Fig. 24, this minimal choice shows a physical 
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Figure 24: A plot of Eq. (|99| ) with Stt /g Nc (along the vertical axis) against u. 

bound state (actually two) with < n < 1 for the coupling greater than some critical value, 
g^Nc/Sir'^ > 18.28. Such a bound state would be significant because it would mean that the short- 
lived magnetic scalar, which we have introduced as a device, can gain longevity at strong coupling, 
so that it can play the role of a spin gluon in larger diagrams. 



28 



There is aiso a soiution for negative u for all couplings. From the interpretation u = we 
see that these solutions would correspond to complex energies with imaginary part zbvrz/a. In fact 
the relation between u and E is fundamentally ambiguous by the imaginary amount 2iTin / a simply 
due to the discretization of r. We presume that this ambiguity is simply a lattice artifact, and it 
seems likely that the solutions with negative u are also artifacts. However, the ultimate proof of 
their artifactual nature must await a more complete understanding of the continuum limit, which 
will involve taking M ^ oo as well as a — > 0. 

5.2 Strong Coupling M = 4 Scalar Fishnet 

Clearly we need to be able to deal with large values of M if we are to understand QCD. In this 
section, to better understand what this will involve, we look more closely at the M = 4 sector of 
the simpler paradigm scalar theory. In particular we would like to explore the effect of the K^-term 
on the l.h.s. of the strong coupling Eq. (|6T|). Recall that we have set 1 — 6i = . We dropped this 
term in the analysis of Sec. ^ because it made the equation intractable for general n. However, for 
the special case of M = 4 (n = M/2 = 2) it is possible to solve this equation. We also note that 
this case was not covered in Sec. ^ (even with the /t^-term removed) because of the limitation to 
odd n. The additional special case of M = 6 will be investigated in the following subsection. 
For M = 4 the strong coupling eigenvalue equation reads 



^2^-(q?+ql)/2To f ^ + - ki , qT/2 + ki - k2)e-('^?+'^i)/^« , (100) 

where = qi + q2 . After the change of variables (k^^ = ki + k2 and k = k2 — ki ) we can integrate 
with respect to k-p on the r.h.s. What remains is 



(^t2_^2^2g-(qf+qi)/2To^^^(q^^q^) ^ 

-{c?+qi)/2To j + k, q^/2 - k)e-'^V2To_ (iq^) 
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If we now perform a variable transformation on the q's (q^ = qi + q2 and q = (qi — q2)/2) 
then the equation becomes 

{t^ _ ^2^2g-(q|+4q^)/4To^ + q, qT/2 - q) = 

A!g-(q2+4q2)/4To /" (qT/2 + k, qT/2 - k)e-'''/2To , (102) 

2 J vrTo 

We can absorb the qy dependence by scaling t = to exp(— q|,/8To), in effect going to the center of 
mass frame. 

For the case of > k?X'^, we can eliminate the ip2 dependence by manipulating this equation 
and integrating with respect to cfq. The result 

is a transcendental equation for the eigenvalue, as is readily seen by direct evaluation of the integral: 

A In 4^-1, (104) 
2k t-K ^ ^ 
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where t = Iq/X and we have assumed k? > 0. If < we write k = ia and the equation becomes 



1 - = - tan"^ J . (105) 
It is immediately clear that a solution exists in this case only for a"^ = —k^ < 1. 



If we analyze (104) we see that by varying i between k and oo the r.h.s. takes on values between 
and oo, thus there is a solution to this equation for any value of k. For the special case of k = 0, 
a more careful analysis of this equation also yields a solution. For a given eigenvalue solution to 



(|104| ) the eigenfunction is given by 

where the integral on the r.h.s. represents a number that can be fixed by normalization. 

If we refer back to (102) we see that if P < k^, then a delta function term restricting momentum 



to energy shell may be added. For s-waves (non s- waves are free), we then have 

where the coefficient of the delta function can be fixed via the following equation which relates A 
to the normalization of the wavefunction. 
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In this case, P < k^, there is no restriction on the value of i, thus this solution corresponds to 
the continuum. To summarize the spectrum of tg includes a discrete s-wave bound state and a 
continuum for tg < k^X^. Clearly the discrete state will not change drastically in the limit k ^ 0, 
but the continuum would be dramatically squeezed to a set of measure zero in this limit. 

5.3 Strong Coupling M = 6 Scalar Fishnet 

The starting point for evaluating the strong coupling M = 6 bit scalar fishnet diagrams is Eq. (^ij) . 
For n = M/2 = 3 the equation to solve is 
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V'3((qi + q2)/2 + k2i, (q2 + q3)/2 + kgs, (q3 + qi)/2 + ki3)e-(l'i+''2+ki)/7^o^ (igg) 

where we have used kjj = kj — kj. By working in the center of mass frame, we can replace 
q3 — > — qi — q2 and the equation becomes 

('t2_^2;^3g-(q?+qi+q,.q.)/To\^3(q^^q^^_q^_q^) ^ ^-{ce^+c^l+^.-c,,) /T, fd^d^fh 
V / J vrTo vrTo ttTq 

V^3(qi/2 + q2/2 + k2i, -qi/2 + kgz, -q2/2 + ki3)e-('^?+^2+ki)/To _ (^q) 

If we change the integration variables k2 and ks on the r.h.s. to 

pi =qi/2 + q2/2 + k2i, p2 = -qi/2 + k32, (111) 
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then we see that Vs in the integrand on the r.h.s. is independent of ki, so the Gaussian integral 
over ki may be trivially performed by completing the square. 
Once this has been done the result is 
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In this equation we see that except for the last term in the exponential on the r.h.s. this integral 
equation only depends on the scalar quantity, qf + q2 + qi • q2- This is not surprising since this 
scalar quantity is proportional to qi + qi + qi and is the only cyclically invariant (qj — > qj+i) scalar 
of order q^. We next perform another change of variables, 
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remembering that q^, pi and now Uj, v, are Euclidean 2- vectors. With these substitutions our 
equation becomes 
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If we rescale U2 and V2 by a factor of -\/3, then we can make the 0(4) symmetry manifest (up to 
the last four terms in the exponential on the r.h.s.) by combining ui and U2 (vi and V2) into a 
Euclidean 4- vector U (V). Thus the equation may be written as 

(i2 _ ^2;^3^-3U^^ ^^^U) = l^^e-^uV2y■^4^^3(v)e-2V^+2V-7^u^ 

= I^g-7UV2 j c^VV'3(7^V)e-2V^+2V.u^ (115) 

where Tl is the real orthogonal 0(4) rotation: 
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We note that = -I. We can search for an 0(4) invariant solution to this equation which is 
a function only of the length, ||U||. Although this will not yield the most general cigenstate, it is 
expected to include the ground state. Plugging in this ansatz, the equation simplifies to 

(t''-K''Xh-^''')MU) = ^e-'^'/^ jv^dViP^{V)e-^^' j d^W-"". (117) 
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where U = ||U|| and V = ||V||. The angular integral on the r.h.s. may be evaluated with standard 
techniques, 



2tt^Ii{2UV) 
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where Ii{x) is the modified Bessel function regular at x = 0. 
Thus the integral equation to solve is 

— 7U^/2 c 

t^ijsiU) = K^X\-^^'iP3iU) + 8\^ '' ^ V^dV hi2UV)My)e-^^' ■ (119) 
The first thing to note is that for = the eigenfunction solution is a Gaussian of the form 

^P^^U) = e~^^\ where g= (120) 
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For this solution the corresponding eigenvalue is 

[11 + VwE) 

Both of these match the values predicted by Equations ( |76| ) and (^) for n = 3. For away from 
zero we can solve this integral equation by means of an iterative procedure. Starting with the 
solution for = we can iterate the r.h.s. of ( |119| ) repeatedly. This is a convenient way of solving 
this equation since the functions generated by the integral on the r.h.s. are always Gaussian. Thus 
at each iteration step the solution will be of the form 

MU) = ^Cne-''"''\ (122) 

n 

with the number of terms in the sum doubling after each iteration. For a solution of this form it 
can be shown that the corresponding eigenvalue, t^, is given by 

,.^^. + iE^W(2W_ (123) 

We have tested this iteration procedure numerically and for values of small (k^ < 0.3) we see 
that the wavefunction (eigenvalue) converges to a well-defined function (value). For larger (closer 
to 1) this becomes murkier as one needs a lot more iterations for the convergence to a value distinct 
from to be evident. Another interesting phenomenon is that for < —0.11, then becomes 
negative indicating that no physical solution exists. 



6 Conclusion 

In this article, we have refined and extended an approach, proposed in the late seventies, to obtain 
the large Nc limit of QCD by directly summing the planar diagrams which survive. The basic tool 
is to define the planar diagrams using light-front space-time coordinates for which ix^ and the 
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carried by each gluon are discretized. This effectively digitizes the sum of diagrams, a first step 
toward a numerical evaluation. It also regulates the usual divergences of Feynman diagrams. We 
identified several shortcomings of the discretized model of QCD attempted in Q, and we proposed 
an improved formulation which at least mitigates, and might well overcome, these defects. 

Discretization enables a formal strong 't Hooft coupling limit of the sum of diagrams. A major 
disadvantage of the discretization of ||5| was that this formal limit suppressed the cubic gluonic 
interaction essential for the "anti-ferromagnetic" ordering of glueball mass levels: the dominant 
quartic interaction ordered levels ferromagnetically. Our new discretized model replaces the quartic 
interactions by the exchange of two kinds of fictitious "short-lived" scalars, so that all interactions 
can compete on an equal footing in the strong coupling limit. The ambiguities inherent in such a 
replacement can also be exploited to remove unwanted symmetry violations induced by the usual 
ultraviolet divergences present in the continuum limit. 

Having defined our discretized model, we explored its physical properties in several ways. We 
first studied the nature of weak coupling perturbation theory by calculating the gluon self energy 
to one loop order, regaining the known continuum answer. This calculation showed how the dis- 
cretization regulates ultraviolet divergences, and how the ambiguities in the model begin to be fixed 
by the restoration of Poincare invariance. Although we have not done a two loop calculation, there 
is sufficient flexibility in these ambiguities to hope to achieve Poincare invariance to all orders in 
perturbation theory. The discretized model can also be studied in the strong coupling limit, but 
in this article we just began this study for QCD by looking only at states with very small total 

= Mm for M = 2, where the dynamics is so drastically simplified that it can be solved exactly. 
We defer to a future publication studies of QCD at M = 3 and higher. The continuum limit, of 
course, will require M — > oo. 

As a warmup for going to larger values of M, we evaluated the strong coupling limit in a 
paradigm matrix scalar field theory with only cubic interactions. Not surprisingly, the bosonic 
light-cone string was obtained. Although this paradigm model yielded some useful insights into 
the nature of large planar diagrams, we stressed that the corresponding QCD calculation will have 
profound differences: for one thing the gluons carry spin, and for another their interactions show 
both repulsion and attraction depending on the quantum numbers of the channel. In contrast, the 
interactions of the scalar theory are exclusively attractive. Because of this, the strong coupling 
limit forced the p"*" carried by each scalar quantum to be minimal, i.e. one discretized unit m. This 
circumstance prevented a "Liouville" degree of freedom, associated with collective fluctuations of 
the P"*" distribution among the scalar quanta, from arising. Thus the limit must be interpreted as 
a critical string theory. 

The diversity of interaction signs of QCD will obviously complicate this outcome. It is possible, 
and a major focus for future study, that cancellations deemphasize the contributions where all 
quanta carry the minimum to such a degree that a collective Liouville fleld emerges. Then the 
strong coupling limit might be a subcritical version of one of the existing string models. If so, 
the Liouville world-sheet field could be thought of as a fifth dimension, and the dual description 
of our model as a field theory at weak coupling and a subcritical string theory at strong coupling 
would resemble the anti de Sitter gravity (AdS) / conformal field theory (CFT) duality of [§,|l7, 18|. 
Another logical possibility, though, is that the strong coupling limit of large Nc QCD is actually 
a novel critical string theory with critical dimension 4. Of course, it could also turn out that the 
attempt to reach a reasonable Poincare invariant strong coupling limit of large Nc QCD simply 
fails. After all, continuum QCD is, strictly speaking, not an infinite coupling theory in any sense 
of the word. The coupling is scale dependent and corresponds to no tunable parameter at all. The 
strong coupling limit, as everyone knows, describes the discretized model and can vary wildly from 
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one discretization to another. 

Much has been said about the "holographic" nature of the duahty mentioned above. We would 
like to conclude with a few comments about this. The hologram metaphor was invented by 't 
Hooft to describe a possible resolution of the "information loss paradox" of quantum black 
holes. Since the horizon of a black hole is two dimensional, it should be possible to describe all of 
three dimensional physics by a two dimensional quantum theory. The discretized model we have 
presented here is not holographic in this sense. The transverse space of a light front is indeed 
two dimensional, but the third longitudinal dimension has not been eliminated: it is present in the 
disguised form of a variable Newtonian mass Mm for each gluon. However, the model is holographic 
in the higher dimensional sense described by Witten The "fundamental" discretized model is 
3 + 1 dimensional, 2 transverse dimensions, variable and rr"*". However in the strong coupling 
limit we expect 4 + 1 dimensions: the x~ of light-cone string should emerge as a function of the 
transverse and Liouville degrees of freedom. Holography in 't Hooft 's sense would require a more 
profound circumstance: there should be no Liouville field and the variable of each gluon must 
itself be a mere collective effect. For example, the gluon with M units of might be thought of 
as a bound system of M minimal p'^ "bits" (2^. In that case, the model presented here would just 
be a stepping stone toward that more fundamental theory. 
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A Appendix: Alternate Discretization 

In this appendix we explore the ramifications of the alternate discretization Eq. (|8|) of D . The 
bare propagator in energy representation is 

The self-energy parts n'^',H'^"^ and H^^ will of course have different values in this discretization, 
but the decompositions (|l^ remain valid. Under the assumption that H'/ = U'^ = Hi, the relations 
of the exact propagators to the H's are identical to (|15|) except for D for which the relation is 

^ To r n(e"^ - 1 - Q72Mro) u 

M2 [l - n - TouUi{e''^ - 1 - Q'^/2MTo) 2MTo 1-u- uU2i/2M 

The only parts of the one loop self energy calculation affected by the different discretization 
are the two diagrams. Fig. |T2|, which have a D propagator as one of the internal lines. The 
evaluation is quite different for this discretization because the completion of squares in the second 
term of Eq. (^) leads to different factors than the first term. The contribution of the first term 
involves the exponent 

k fp' , {Cl-P? \_ k f M{p-lQ/M)^ Q2\ 
2To \ I M-l J 2To V l{M-l) M J ' ^ ' ^ 

whereas the contribution of the second term involves 

(k - i)p2 ^ fc(Q-p)' = — f + ^ (p - _J^Q__\ ' + k{k -l)Q^ \ 



2TqI 2Tq{M-1) 2Tq\ l{M-l) V M{k-l)+lJ M{k-l) + l 
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for one of the two diagrams, and for the other a similar expression with p ^ Q — p and I 
Thus the two contribute equally to IIj^j: 



for k > 1 and by 



47r2 



(2M - 0" ^_fcQ2/2AfTo 



„fc«QV2MTo(Af(A:-l)+/) 

M(A;-1) + / 



1 



MA: 



M-L 



(A.5) 



g^N, (2M - l)^ e-QV2MTo 
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M 



(A.6) 



for k = 1. Translating to energy representation gives the more compact 



frAV 
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' Mk + l 
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(A.7) 



Adding the result of the unchanged first thirteen diagrams to this and subtracting Q2ni/2 gives 
112 for this discretization: 



g'^Nc 13M2 - 12M 
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M M^J To ^ ' 
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k=l 



487r2 

g/{A:+l)Q2/2MTo(Mfc+0 

Mk+l Hk 



(Ai 



The violation of Galilei invariance caused by this alternate discretization is apparent from the non- 
polynomial dependence on Q2. However, one can easily see that each power of comes with an 
accompanying power of 1/M. At most one power of M is supplied by the prefactors, so all powers 
of Q2 higher than the first are irrelevant in the continuum limit. 

In order to compare our results to the continuum calculation of Ref. ||8|, we must examine 
the limits n — > 1 and M — > oo. The behavior of the first two terms of Eq. ( A.8| ) in this limit is 
transparent once we use the identity (|40|). The continuum limit of the last term requires a bit 



more analysis. First, as mentioned in the previous paragraph, we only need keep two terms in the 
expansion of the exponential: 
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The sums over / can be approximated using the Euler-Maclaurin summation formula, see Eq. (^3|), 
as long as F is not singular at the endpoints of integration. Clearly terms with 1 // in the summand 
must be treated separately, which is easily handled using (35). Applying these formulae, we find 
(for large M but arbitrary u). 
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The continuum limit also requires 1 — u ~ Q"^ /2MTq = (Q^ — 2inME)/2MTQ, so u may be set to 
unity in all nonsingular terms without a prefactor of M. Then the above terms simplify to 
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where we have defined 
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which can be numerical evaluated: 

a ^ 1.188, /3 « 1.991, (5 « 0.633. 
Putting all this together, we get for the continuum limit. 
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Remembering that our 112 is a factor of —Q'^/Tq times that defined in Ref. [Q, we find agreement 
for the coefficient of In , provided we identify 2MTo with A2 . The first two groups of terms on the 
r.h.s. of Eq. ( [A.16| ) violate important symmetries and must be removed by explicit counter-terms. 
The term linear in Q2 violates Galilei invariance and the momentum independent terms imply a 
finite gluon mass squared in perturbation theory, thus violating Poincare invariance. Note that 
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the necessary counter terms are low order polynomials in both the transverse momentum and in 
the discretized P"*". Without the tadpole contributions the discretization used in the text would 
have required counter-terms with logarithmic M dependence. But we found (at least at one-loop) 
that the tadpoles could be designed to eliminate the need for counter-terms. Then the fact that it 
preserves Galilei invariance makes it the superior choice. 
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